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Recent advances in experimental techniques now permit to measure the Casimir force with un- 
precedented precision. In order to achieve a comparable precision in the theoretical prediction of 
the force, it is necessary to accurately determine the electric permittivity of the materials consti- 
tuting the plates along the imaginary frequency axis. The latter quantity is not directly accessible 
to experiments, but it can be determined via dispersion relations from experimental optical data. 
In the experimentally important case of conductors, however, a serious drawback of the standard 
dispersion relations commonly used for this purpose, is their strong dependence on the chosen low- 
frequency extrapolation of the experimental optical data, which introduces a significant and not 
easily controllable uncertainty in the result. In this paper we show that a simple modification of the 
standard dispersion relations, involving suitable analytic window functions, resolves this difficulty, 
making it possible to reliably determine the electric permittivity at imaginary frequencies solely 
using experimental optical data in the frequency interval where they are available, without any need 
of uncontrolled data extrapolations. 

PACS numbers: 05.30.-d, 77.22. Ch, 12.20.Ds 
Keywords: Casimir, dispersion relations. 
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I. INTRODUCTION 

One of the most intriguing predictions of Quantum 
electrodynamics is the existence of irreducible vacuum 
fluctuations of the electromagnetic (e.m.) field. It was 
Casimir's fundamental discovery [l| to realize that this 
purely quantum phenomenon was not confined to the 
atomic scale, as in the Lamb shift, but would rather man- 
ifest itself also at the macroscopic scale, in the form of 
a force of attraction between two discharged plates. For 
the idealized case of two perfectly reflecting plane-parallel 
plates at zero temperature, placed at a distance a in vac- 
uum, Casimir obtained the following remarkably simple 
estimate of the force per unit area 



F a = 



ir 2 hc 
240 a 4 



(1) 



An important step forward was made a few years later 
by Lifshitz and co-workers Q, who obtained a formula 
for the force between two homogeneous dielectric plane- 
parallel slabs, at finite temperature. In this theory, of 
macroscopic character, the material properties of the 
slabs were fully characterized in terms of the respec- 
tive frequency dependent electric permittivities e(w), ac- 
counting for the dispersive and dissipative properties of 
real materials. In this way, it was possible for the first 
time to investigate the influence of material properties 
on the magnitude of the Casimir force. 

Over ten years ago, a series of brilliant experiments 
0, 0] exploiting modern experimental techniques pro- 
vided the definitive demonstration of the Casimir effect. 
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These now historical experiments spurred enormous in- 
terest in the Casimir effect, and were soon followed by 
many other experiments. The subsequent experiments 
were aimed at diverse objectives. Some of them explored 
new geometries: while the works 0, H| used a sphere- 
plate setup, the original planar geometry investigated by 
Casimir was adopted in the experiment [9 , and a setup 
with crossed cylinders was considered in [6|. The impor- 
tant issue of the non trivial geometry dependence of the 
Casimir effect is also being pursued experimentally, us- 
ing elaborate micro- patterned surfaces ■ Other exper- 
iments aimed at demonstrating new possible uses of the 
Casimir force, like for example the actuation of microma- 
chines 0, or at demonstrating the possibility of a large 
modulation of the Casimir force [H, [l(| , which could also 
result in interesting technological applications. There are 
also experiments using superconducting Casimir cavities, 
that aim at measuring the change of the Casimir energy 
across the superconducting phase transition [ll[. The 
experiments performed in the last ten years are just too 
numerous to mention them all here. For an updated ac- 
count we refer the reader to the very recent review paper 

Apart from exploring new manifestations of the 
Casimir effect, a large experimental effort is presently be- 
ing made also to increase the precision of Casimir force 
measurements, in simple geometries. Already in the early 
experiment Q a precision upto one percent was obtained. 
More recently, a series of experiments with microtorsional 
oscillators 1 1 31 ] reached an amazing precision of 0.2 per- 
cent. The reader may wonder what is the interest in 
achieving such a high precision in this kind of experi- 
ments. There are several reasons why this is important. 
On one hand, in the theory of dispersion forces puzzling 
conceptual problems have recently emerged that are con- 
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nected with the contribution of free charges to the ther- 
mal Casimir force, whose resolution crucially depends on 
the precision of the theory-experiment comparison fl2l ]. 
On the other hand, the ability to accurately determine 
the Casimir force is also important for the purpose of ob- 
taining stronger constraints on hypothetical long-range 
forces predicted by certain theoretical scenarios going be- 
yond the Standard Model of particle physics [12j. 

The remarkable precision achieved in the most recent 
experiments poses a challenging demand on the theorist: 
is it possible to predict the magnitude of the Casimir 
force with a comparable level of precision, say of one 
percent? Assessing the theoretical error affecting present 
estimates of the Casimir force is a difficult problem in- 
deed, because many different factors must be taken into 
account [l2[ . Consider the typical experimental setting of 
most of the current experiments, where the Casimir force 
is measured between two bodies covered with gold, placed 
in vacuum at a distance of a (few) hundred nanometers. 
In this separation range, the main factor to consider is 
the finite penetration depth of electromagnetic fields into 
the gold layer [25| . resulting from the finite conductivity 
of gold. The tool to analyze the influence of such mate- 
rial properties as the conductivity on the Casimir effect 
is provided by Lifshitz theory [2j. This theory shows 
that for a separation of 100 nm, the finite conductivity 
of gold determines a reduction in the magnitude of the 
Casimir force of about fifty percent in comparison with 
the perfect metal result [14]. Much smaller corrections, 
that must nevertheless be considered if the force is to 
be estimated with percent precision, arise from the finite 
temperature of the plates and from their surface rough- 
ness. Moreover, geometric effects resulting from the ac- 
tual shape of the plates should be considered. We should 
also mention that the magnitude of residual electrostatic 
forces between the plates, resulting from contact poten- 
tials and patch effects, must be carefully accounted for. 
For a discussion of all these issues, which received much 
attention in the recent literature on the Casimir effect, 
we again address the reader to Ref. ij.2j . See also the 
recent work [15j . 

In this paper, we focus our attention on the influence 
of the optical properties of the plates which, as explained 
above, is by far the most relevant factor to consider. As 
we pointed out earlier, in Lifshitz theory the optical prop- 
erties of the plates enter via the frequency-dependent 
electric permittivity e(uj) of the material constituting 
plates. In order to obtain an accurate prediction of the 
force, it is therefore of the utmost importance to use ac- 
curate data for the electric permittivity. The common 
practice adopted in all recent Casimir experiments with 
gold surfaces is to use tabulated data for gold (most of the 
times those quoted in Refs. 16]), suitably extrapolated 
at low frequencies, where optical data are not available, 
by simple analytic models (like the Drude model or the 
so-called generalized plasma model). However, already 
ten years ago Lamoreaux observed [17j that using tabu- 
lated data to obtain an accurate prediction of the Casimir 



force may not be a reliable practice, since optical prop- 
erties of gold films may vary significantly from sample to 
sample, depending on the conditions of deposition. The 
same author stressed the importance of measuring the 
optical data of the films actually used in the force mea- 
surements, in the frequency range that is relevant for the 
Casimir force. The importance of this point was further 
stressed in [l8j and received clear experimental support 
in a recent paper pjj], where the optical properties of 
several gold films of different thicknesses, and prepared 
by different procedures, were measured ellipsometrically 
in a wide range of wavelengths, from 0.14 to 33 microns, 
and it was found that the frequency dependent electric 
permittivity changes significantly from sample to sam- 
ple. By using the zero-temperature Lifshitz formula, the 
authors estimated that the observed sample dependence 
of the electric permittivity implies a variation in the the- 
oretical value of the Casimir force, from one sample to 
another, easily as large as ten percent, for separations 
around 100 nm. It was concluded that in order to achieve 
a theoretical accuracy better than ten percent in the pre- 
diction of the Casimir force, it is necessary to determine 
the optical properties of the films actually used in the 
experiment of interest. 

The aim of this paper is to improve the mathemati- 
cal procedure that is actually needed to obtain reliable 
estimates of the Casimir force, starting from experimen- 
tal optical data on the material of the plates, like those 
presented in Ref. [l9j. The necessity of such an im- 
provement stems from the very simple and unavoidable 
fact that experimental optical data are never available 
in the entire frequency domain, but are always restricted 
to a finite frequency interval u> m i n < u> < u; max . To see 
why this constitutes problem we recall that Lifshitz for- 
mula, routinely used to interpret current experiments, 
expresses the Casimir force between two parallel plates 
as an integral over imaginary frequencies i£ of a quan- 
tity involving the dielectric permittivities of the plates 
e(i£). For finite temperature, the continuous frequency 
integration is replaced by a sum over discrete so-called 
Matsubara frequencies ui n = i£„, where £ n = 2TmltBT/h, 
with n a non- negative integer, and T the temperature of 
the plates. In any case, whatever the temperature, one 
needs to evaluate the permittivity of the plates for cer- 
tain imaginary frequencies. We note that, in principle, 
recourse to imaginary frequencies is not mandatory be- 
cause it is possible to rewrite Lifshitz formula in a math- 
ematically equivalent form, involving an integral over the 
real frequency axis. In this case however the integrand 
becomes a rapidly oscillating function of the frequency, 
which hampers any possibility of numerical evaluation. 
In practice, the real-frequency form of Lifshitz formula is 
never used, and only its imaginary-frequency version is 
considered. We remark that occurrence of imaginary fre- 
quencies in the expression of the Casimir force, is a gen- 
eral feature of all recent formalisms, extending Lifshitz 
theory to non-planar geometries |20l - t22j | . The problem is 
that the electric permittivity e(i£) at imaginary frequen- 
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cies cannot be measured directly by any experiment. The 
only way to determine it by means of dispersion relations, 
which allow to express e(i£) in terms of the observable 
real- frequency electric permittivity e(w). In the standard 
version of dispersion relations [2|, adopted so far in all 
works on the Casimir effect, e(i£) — 1 is expressed in terms 
of an integral of a quantity involving the imaginary part 
e"(w) of the electric permittivity: 

The above formula shows that, in principle, a determina- 
tion of e(i£) requires knowledge of e"(w) at all frequencies 
while, as we said earlier, optical data are available only in 
some interval w m i n < to < w max . In practice, the problem 
is not so serious on the high frequency side, because the 
fall-off properties of e"(u) at high frequencies, together 
with the ui 2 factor in the denominator of the integrand, 
ensure that the error made by truncating the integral at 
a suitably large frequency oj max is small, provided that 
large enough. Typically, an w max larger than, say, 
15c/ (2a), is good enough for practical purposes. Things 
are not so easy though on the low frequency side. In 
the case of insulators, optical data are typically available 
until frequencies uj m - lu much smaller than the frequen- 
cies of all resonances of the medium. Because of this, 
e"(w) is almost zero for uj < uj lri \ n , and therefore the 
error made by truncating the integral at w m i n is again 
negligible. Problems arise however in the case of ohmic 
conductors, because then e"(uj) has al/u singularity at 
uj = 0. As a result e"(w) becomes extremely large at low 
frequencies, in such a way that the integral in Eq. ([2]) 
receives a very large contribution from low frequencies. 
For typical values of w m i n that can be reached in practice 
(for example for gold, the tabulated data in [16| begin 
at cj m i n = 125 meV/h, while the data of start at 38 
meV/?i) truncation of the integral at w m i n results in a 
large error. The traditional remedy to this problem is to 
make some analytical extrapolation of the data, typically 
based on Drude model fits of the low-frequency region of 
data, from uj m i n to zero, and then use the extrapolation 
to estimate the contribution of the integral in the interval 
< uj < uj m - m where data are not directly available. It 
is important to observe that this contribution is us ually 
very large. For example, even in the case of Ref. (19| . 
the relative contribution of the extrapolation is about 
fifty percent of the total value of the integral, in the en- 
tire range of imaginary frequencies that are needed for 
estimating the Casimir force. 

Clearly, this procedure is not very satisfying. The use 
of analytical extrapolations of the data introduces an un- 
certainty in the obtained values of e(i£), that is not easy 
to quantify. The result may in fact depend a lot on the 
form of the extrapolation, and there is no guarantee that 
the chosen expression is good enough. Consider for ex- 
ample Ref. uM, which constitutes the most accurate 
existing work on this problem. It was found there that 
the simple Drude model does not fit so well the data 



of all samples, making it necessary to improve it by the 
inclusion of an additional Lorentz oscillator. Moreover, 
it was found that for each sample the Drude parame- 
ters extracted from the data depended on the used fit- 
ting procedure, and were inconsistent which each other 
within the estimated errors, which is again an indication 
of the probable inadequacy of the analytical expression 
chosen for the interpolation. This state of things led us 
to investigate if it possible to determine accurately e(i£) 
solely on the basis of available optical data, without mak- 
ing recourse to data extrapolations. We shall see below 
that this is indeed possible, provided that Eq. ^ is 
suitably modified, in a way that involves multiplying the 
integrand by an appropriate analytical window function 
f(oj), which suppresses the contribution of frequencies 
not belonging to the interval w m i n < uj < ui max . As a 
result of this modification, the error made by truncating 
the integral to the frequency range uj m - m < uj < 0J max 
can be made negligible at both ends of the integration 
domain, rendering unnecessary any extrapolation of the 
optical data outside the interval where they are avail- 
able. The procedure outlined in this paper should allow 
to better evaluate the theoretical uncertainty of Casimir 
force estimates resulting from experimental errors in the 
optical data. 

The plan of the paper is as follows: in Sec. II we de- 
rive a generalized dispersion relation for e(i£), involving 
analytic window functions f(z), and we provide a simple 
choice for the window functions. In Sec III we present 
the results of a numerical simulation of our window func- 
tions, for the experimentally relevant case of gold, and in 
Sec IV we estimate numerically the error on the Casimir 
pressure resulting from the use of our window functions. 
Sec V contains our conclusions and a discussion of the 
results. 



II. GENERALIZED DISPERSION RELATIONS 
WITH WINDOW-FUNCTIONS 

As it it well known Q, analyticity properties satisfied 
by the electric permittivity e(w) of any causal medium 
(and more in general by any causal response function, the 
magnetic permeability fJ>(uj) being another example) im- 
ply certain integral relations between the real part e'(u>) 
and imaginary part e"(ui) of e(w), known as Kramers- 
Kronig or dispersion relations. The dispersion relation 
of interest to us is the one that permits to express the 
value e(i£) of the response function at some imaginary 
frequency i£ in terms of an integral along the positive fre- 
quency axis, involving e"(u;). It is convenient to briefly 
review here the simple derivation of this important re- 
sult, which is an easy exercise in contour integration. 
For our purposes, it is more convenient to start from 
an arbitrary complex function u(z), with the following 
properties: u(z) is analytic in the upper complex plane 
C + = {z : Im(z) > 0}, fall's off to zero for large \z\ 
like some power of \z\, and admits at most a simple pole 
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at 10 = 0. Consider now the closed integration contour 
r obtained by closing in the upper complex plane the 
positively oriented real axis, and let zq be any complex 
number in C + . It is then a simple matter to verify the 
identity: 



dz 



z u(z) 



mu(zo) 



(3) 



The assumed fall-off property of u(z) ensures that the 
half-circle of infinite radius forming T contributes nothing 
to the integral, and then from Eq. ([3]) we find: 



u(z ) = T~ 



°° , LUu(ui) 
dUJ^: 7T 



(4) 



Consider now a purely imaginary complex number zq — 
i£, and assume in addition that along the real axis u(to) 
satisfies the symmetry property u(— ui) = u*(u>). From 
Eq. Q we then find: 



00 , wu"(w) 
dui- 



uj- 



(5) 



which is the desired result. 

The standard dispersion relation Eq. ([2]) used to com- 
pute the electric permittivity for imaginary frequencies 
is a special case of the above relation, corresponding to 
choosing u(z) = e(z) — 1. We note that Eq. ([2]) is valid 
both for insulators, which have a finite permittivity at 
zero frequency, as well as for ohmic conductors, whose 
permittivity has al/w singularity in the origin. As we 
explained in the introduction Eq. ([2]), even though per- 
fectly correct from a mathematical standpoint, has seri- 
ous drawbacks, when it is used to numerically estimate 
e(i£) for ohmic conductors, starting from optical data 
available only in some interval w m i n < u> < w max , be- 
cause the integral on the r.h.s. of Eq. © receives a large 
contribution from frequencies near zero, where data are 
not available. This difficulty can however be overcome 
in a very simple way, as we now explain. Consider a 
window function f(z), enjoying the following properties: 
f(z) is analytic in C + , it has no poles in C + except possi- 
bly a simple pole at infinity, and satisfies the symmetry 
property 



/(-**) = /*(*) 



(6) 



Consider now Eq. for u(z) = f(z)(e(z) — 1). Since for 
any medium (e(z) — 1) falls off like z~ 2 at infinity [23j . 
the quantity u(z) falls off at least like z _1 at infinity, 
and it satisfies all the properties required for Eq. (|5|) to 
hold. For any £ such that f(i£) ^ 0, we then obtain the 
following generalized dispersion relation: 



em - 1 = 



dui- 



Im[/H(eH - 1)] 



(7) 

We note that the above relation constitutes an exact re- 
sult, generalizing the standard dispersion relation Eq. 



([2]), to which it reduces with the choice f(z) = 1. An- 
other form of dispersion relation, frequently used in the 
case of conductors or superconductors [H|,[2J] is obtained 
by taking f(z) = iz into Eq. ([7]). Recalling the relation 
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u> 



it reads: 



du- 



rim [a(tu) 



(8) 



(9) 



The above form is especially convenient in the case of 
superconductors, because it avoids the 5(v) singularity 
characteri zing the real part of the conductivity of these 
materials [24j |. 

We observe now, and this is the key point, that there 
is no reason to restrict the choice of the function f(z) 
to these two possibilities. Indeed, we can take advan- 
tage of the freedom in the choice of f(z), to suppress 
the unwanted contribution of low frequencies (as well as 
of high frequencies), where experimental data on e(o>) 
are not available. In order to do that, it is sufficient to 
choose a window function that goes to zero fast enough 
for lj — > 0, as well as for uj — ¥ oo. A convenient family of 
window functions which do the job is the following: 



f{z) = Az^ +1 



1 



f 



(z-w) 2 i +1 (z + w*) 2 ^ 1 



(10) 



where w is an arbitrary complex number such that 
Im(w) < 0, and p and q are integers such that p < q. 
The constant A is an irrelevant arbitrary normalization 
constant, that drops out from the generalized dispersion 
formula Eq. ([7]). As we see, in the limit z — > 0, these 
functions vanish like z 2p+1 , and therefore by taking suf- 
ficiently large values for p we can obtain suppression of 
low frequencies to any desired level. On the other hand, 
for z — > oo, f(z) vanishes like z 2 ^ p ~ q \ and therefore 
by taking sufficiently large values of q, we can obtain 
suppression of high frequencies. Moreover, by suitably 
choosing the free parameter w, we can also adjust the 
range of frequencies that effectively contribute to the in- 
tegral on the r.h.s. of Eq. Q. In Figs. 1 and 2 we 
plot the real and imaginary parts (in arbitrary units) of 
our window functions /(w), versus the frequency uj (ex- 
pressed in eV). The two curves displayed correspond to 
the choices p = 1, q = 1 (dashed line) and p = 1, q = 3 
(solid line). In both cases, the parameter w has the value 
w = (1 — 2i) eV/fr. We observe that along the real fre- 
quency axis, our window functions have non- vanishing 
real and imaginary parts. This is not a feature of our 
particular choice of the window functions, but it is an 
unavoidable consequence of our demand of analyticity on 
f(z). Indeed, for real frequencies uj the real and imagi- 
nary parts of f(u>) are related to each other by the usual 
Kramers-Kronig relations Q that hold for the boundary 
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FIG. 1: Real part f'(cj) (in arbitrary units) of the window 
functions in Eq. (I10|l versus frequency to (in eV/h). The 
window parameters are p = 1, q = 2 (dashed line), p = 1, q — 
3 (solid line). In both cases, the parameter w has the value 
w = (1 - 2i) eV/n. 
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FIG. 2: Imaginary part /"(w) (in arbitrary units) of the win- 
dow functions in Eq. (|10[) versus frequency lj (in eV/h). 
The window parameters are p = 1, <j = 2 (dashed line), 
p = 1, q = 3 (solid line). In both cases, the parameter u> 
has the value w = (1 — 2i) eV/h. 
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FIG. 3: Plots (in arbitrary units) of the window functions 
/(i£) in Eq. (JTDJ versus the imaginary frequency £ (in eV/h). 
The window parameters are p = 1, g = 2 (dashed line), p = 
1, g = 3 (solid line). In both cases, the parameter w has the 
io= (1 -2i) eV/ft. 



values of analytic functions. In the case when f[z) van- 
ishes at infinity, they read: 



/'M 



l. 



no 



i r c 

/"(c) = -ip / 



■ w 



(11) 



(12) 



where the symbol P in front of the integrals denotes the 
principal value. These relation show that vanishing of 
f'(uj) implies that of f"(ui) and viceversa, and therefore 
neither f'(tu) nor f"{ui) can be identically zero. By virtue 
of this property of the window functions, it follows from 
Eq. ([7]) that both the real and imaginary parts of e(w) 
are needed to evaluate e(i£) (unless the standard choices 
f(z) = 1 or f(z) = \ z are made). We also note (see Fig. 
1 and 2) that the real and imaginary parts of f(co) do 
not have a definite sign. This feature also is a general 
consequence of our key demand that f(z) vanishes in the 
origin, as it can be seen by taking u = in Eqs. (fTTj) and 
(1T21) . Since the l.h.s. of both equations are required to 
vanish, the integrand on the r.h.s. cannot have a definite 
sign. Finally, in Fig 3 we show plots of two of our win- 
dow functions versus the imaginary frequency £, 
expressed in eV, for the same two choices of parameters 
of Fig. 1 and 2. It is important to observe that the win- 
dow functions f(z) are real along the imaginary axis (as 
it must be, as a consequence of the symmetry property 
Eq. (j6|)). However, the sign of is not definite, and 
as a result of this admits zeros along the imaginary 
axis. When using Eq. {7} for estimating e(i£) it is then 
important to choose the window function such that none 
of its zeroes coincides with the value of £ for which e(i£) 
is being estimated. 



III. A NUMERICAL SIMULATION 

In this Section, we perform a simple simulation to test 
the degree of accuracy with which the quantity e(i£) can 
be reconstructed using our window functions, starting 
from data on e(w) referring to a finite frequency interval. 
To do that we can proceed as follows. 

According to the standard dispersion relation Eq. © , 
the quantity e(i£) — 1 is equal to the integral on the r.h.s. 
of Eq. (J5J. Following Refs. QJ, [Hj], we can split this 
integral into three pieces, as follows: 



2 r , uie"(ui) 
duj— — 



71" Jo w 
where we set: 



^Low(0 — — 



/low(£)+W£)+/high(£) , (13) 

u)e"(oj) 



2 

7T ./,, 



duj- 



wo 



did 



me"(uj) 
'oj 2 + e 



(14) 



(15) 
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and 



T _ 2 Z" 00 J ue«(w) 
4igh(0 = - / duj- 



(16) 



'w 2 + £ 2 ' 
By construction, we obviously have: 

e (i0 - 1 = /iow(0 + WO + W(0 • ( 17 ) 

An analogous split can be performed in the integral on 
the r.h.s. of the other standard dispersion relation in- 
volving the conductivity Eq. (0): 



duj- 



(18) 

with an obvious meaning of the symbols. Again, we have 
the identity: 

e(i£) - 1 = X low (0 + K cxp (0 + K high (0 . (19) 

On the other hand, according to our generalized disper- 
sion relation Eq. 0, the quantity e(i£) — 1 is also equal 
to the integral on the r.h.s. of Eq. (J7J. We can split this 
integral too in a way analogous to Eq. fj 13[) : 



diu- 



Im[/M(e(w) - 1)] 



where we set: 



(21) 

2 /-""■ 



J ixp ?) (0 = / ^-^ 7I Im[/( W )(e(c) - 1)] 



c^ 2 +£ 2 



and 

"'high \<s) ~ 



(22) 



(23) 



Then by construction we also have: 

" 1 = 4^ (0 + ^ (0 + J&g (0 • (24) 

The quantities 2exp(£)j -^exp(C) and Jex P \0 evidently 
represent the contribution of the experimental data. On 

the contrary the quantities Ii ow (£), i^low(£) and 4^(0 
can be determined only by extrapolating the data in the 
low frequency region < lo < w m i a , while determination 

of the quantities 4igh(£). A 4ii g h(£) and (0 is onl y 
possible after we extrapolate the data in the high fre- 
quency interval w max < w < oo. Ideally, we would like 

to have J low (£), WO. Kiow(0 ) *W0> J iow 9) (0 and 
^Mgh (0 as small as possible. 



To see how things work, we can perform a simple sim- 
ulation of real experimental data. We imagine that the 
electric permittivity of gold is described by the follow- 
ing six-oscillators approximation 12], which is known to 
provide a rather good description of the permittivity of 
gold for the frequencies that are relevant to the Casimir 
effect: 



uj(uj + 17) 



3=1 3 



9j 



(25) 



Here, u> p is the plasma frequency and 7 is the relaxation 
frequency for conduction electrons, while the oscillator 
terms describe core electrons. The values of the param- 
eters <7j, ujj and 7j can be found in the second of Refs. 
13]. For ujp and 7 we use the reference values for crys- 
talline bulk samples, ui p = 9 eV/H and 7 = 0.035 eV/Ti. 
Of course with such a simple model for the permittivity 
of gold, there is no need to use dispersion relations to ob- 
tain the expression of e(i£), for this can be simply done 
by the substitution uj — > i£ in the r.h.s. of Eq. (|2"5)) : 



e(i0 



£(e+7) 



3=1 j 



9j 



(26) 



Simulating the real experimental situation, let us pre- 
tend however that we know that the optical data of 
gold are described by Eq. (|2"5j) only in some interval 
w m i„ < uj < w max , and assuming that we do not want to 
make extrapolations of the data outside the experimen- 
tal interval, let us see how well the quantities / exp (0> 

Kexp(Q and jixp (0 defined earlier reconstruct the ex- 
act value of e(i£) — 1 given by Eq. (|26l) . In our simula- 
tion we took ui m i n = 0.038 eV/h (representing the mini- 
mum frequency value for which data for gold films were 
measured in [191 ]]) while for uj max we choose the value 
w max = 30 eV/h. The chosen value of w max is about 
thirty times the characteristic frequency c/(2a) for a sep- 
aration a = 100 nm. The result of our simulation are 
summarized in Figs. 4 and 5. 

In Fig. 4, we report the relative per cent errors 
5i = 100 [1 - 4x P (£n)/(e(i£n) ~ 1)] (black squares) and 
8 K = 100 [1 - K exp (£n)/(e(i£n) - 1)] (grey triangles) 
which are made if the quantities / exp (Cn) or K cxp (^ n ) 
are used, respectively, as estimators of e(i£ n ) — 1. The in- 
teger number on the abscissa labels the Matsubara mode 
f„ = 2unk B T/h (T = 300 K). Only the first sixty modes 
are displayed, which are sufficient to estimate the Casimir 
force at room temperature, for separations larger than 
100 nm, with a precision better than one part in ten 
thousand. As we see, both I eX p(£n) and K exp (t; n ) pro- 
vide a poor approximation to e(i£ n )— 1, with I exp (£,n) per- 
forming somehow better at higher imaginary frequencies, 
and K exp (^ n ) doing better at lower imaginary frequen- 
cies. Indeed, / exp (£n) and K exp (£ n ) suffer from opposite 
problems. On one hand the large error affecting I eKp (£ n ) 
arises mostly from neglect of the large low- frequency con- 
tribution I\ ow (£„ ) , and to a much less extent from neglect 
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Matsubara mode Matsubara mode 



FIG. 4: Numerical simulation of the errors (in percent) in 
the estimate of e(i£n) — 1 for gold, resulting from using the 
quantities / e xp(£n) (black squares) and -R"ex P (£n) (grey trian- 
gles) as estimators, in the hypothesis that data are available 
from Wmin = 0.038 eV/h to w max = 30 eV/h. The integer 
on the abscissa labels the Matsubara mode f n = 2-KnksT '/h 
(T = 300 K). 



FIG. 5: Numerical simulation of the error (in percent) in 
the estimate of e(i£n) — 1 for gold, resulting from using the 
quantity J<S P (£ n ) as a estimator, in the hypothesis that data 
are available from w m i n = 0.038 eV/h to o; max = 30 eV/h. 
The integer on the abscissa labels the Matsubara mode £ n = 
2-KnksT ' i 'h (T = 300 K). Grey triangles are for the window 
function having p = 1, q = 2, black squares for p = 1, q = 3. 
In both cases w = (1 — 2i) eV /h. 



of the high frequency contribution /high(Cn) (The magni- 
tude of the high frequency contribution /high(£n) is less 
than two percent of e(i£ n ) — 1 for all n < 60). The situ- 
ation is quite the opposite in the case of K cxp (^ n ). This 
difference is of course due to the opposite limiting behav- 
iors of the imaginary parts of the permittivity e"(w) in 
the limits lo — > 0, and w — > oo, as compared to those of 
the imaginary part of the conductivity a"(uj). Indeed, for 
v — > 0, e"(oj) diverges like while cr"{uj) approaches 
zero like lo. This explains while the low frequency con- 
tribution Iiow(Cra) is much larger than -Ki ow (£«,)• On the 
other hand, in the limit lu — > oo, e"( w ) vanishes like w~ 3 , 
while ct"(lu) vanishes only like cj . This implies that 
large frequencies are much less of a problem for / C xp(^n) 
than for K exp (^ n ). The conclusion to be drawn from 
these considerations is that, if either of the two stan- 
dard forms Eq. ([2]) or Eq. ([9]) of dispersion relations are 
used, in order to obtain a good estimate of e(i^ n ) — 1, one 
is forced to extrapolate somehow the experimental data 
both to frequencies less than io m i n , and larger than w max . 

We can now consider our windowed dispersion rela- 
tion, Eq. ([7]), with our choice of the window functions 
f(z) in Eq. ([TU| . In Fig. 5, we display the relative per 

cent error 8^ = 100 [1 - j£ p q } (£n)/(e(i£n) - 1)] which 
is made if the quantity jixp^^n) is used as an estimator 
of e(i£n) — 1. We considered two choices of parameters 
for our window functions in Eq. (fit)]) , i.e. p = 1, q = 2 
(grey triangles) and p = 1, q = 3 (black squares). In 
both cases, we took for the parameter w the constant 
value w = (1 — 2i) eV/h (See Figs. 1, 2 and 3). It is ap- 
parent from Fig. 5 that both window functions perform 
very well, for all considered Matsubara modes. The error 

(1 2) 

made by using J CX p (£„) is less than one percent, in ab- 
solute value, while the error made by using Jcxp\^ n ) is 
less than 0.25 percent. The jumps displayed by the rela- 



tive errors in Fig. 5 (around n — 6 for the grey dots, and 
n = 14 for the black ones) correspond to the approximate 
positions of the zeroes of the respective window functions 
/(i£) (see Fig. 3). Such jumps can be easily avoided, fur- 
ther reducing at the same time the error, by making a 
different choice of the free parameter w for each value of 
n. We did not do this here for the sake of simplicity. It is 
clear that in concrete cases one is free to choose for each 
value of n, different values of all the parameters p, q and 
w, in such a way that the error is as small as possible. 

IV. SIMULATION OF THE CASIMIR FORCE 

In this Section, we investigate the performance of our 
window functions with respect to the determination of 
the Casimir force. We consider for simplicity the proto- 
typical case of two identical plane-parallel homogeneous 
and isotropic gold plates, placed in vacuum at a distance 
a. As it is well known, the Casimir force per unit area is 
given by the following Lifshitz formula: 

11 n>0 J ^TE.TM^"^"'*^ / 

(27) 

where the plus sign corresponds to an attraction be- 
tween the plates. In this Equation, the prime over the 
n-sum means that the n — term has to taken with 
a weight one half, T is the temperature, k± denotes 
the magnitude of the projection of the wave- vector onto 
the plane of the plates and q n — \/k^_+ £^/c 2 , where 
= 2ttti hsT/h are the Matsubara frequencies. The 
quantities r a (i£„ , k± ) denote the familiar Fresnel reflec- 
tion coefficients of the slabs for a-polarization, evaluated 
at imaginary frequencies i£„. They have the following 
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FIG. 6: Simulation of the error (in percent), versus plate 
separation (in /im) in the estimate of the Casimir force per 
unit area, between two plane-parallel gold plates in vacuum 
at a temperature T = 300K, resulting from nsing i/(>xp ^(£n) 
as an estimator of e(i£ n ) — 1. The window functions are the 
same as in Fig. 5: the dashed line is for p = 1, q = 2 and the 
solid line for p = 1, q — 3. All values of the other parameters 
are same as in Fig. 5. 



expressions: 



rrE(i£n,fcj 



»"TM(i£n,k_L) 



e(i£n) On - K, 



(28) 



(29) 



where k n = \fk\ + e(i£„)£2/ c 2 . 

We have simulated the error made in the estimate of 
P(a, T) if the estimate of e(i£„) provided by the window- 
approximations Jexp\(,n) is used: 



e(i£„)^l + J&«>te0 



(30) 



again assuming the simple six-oscillator model of Eq. 
(1251) for e(w). The results are summarized in Fig 6, where 

we plot the relative error Sp' q ^ in percent, as a function 
of the separation a (in microns). The window functions 
that have been used are the same as in Fig. 5. We see 
from the figure that already with this simple and not- 
optimized choice of window functions, the error is much 
less than one part in a thousand in the entire range of 
separations considered, from 100 nm to one micron. 



V. CONCLUSIONS AND DISCUSSION 

In recent years, a lot of efforts have been made to 
measure accurately the Casimir force. At the moment 
of this writing, the most precise experiments using gold- 
coated micromechanical oscillators claim a precision bet- 
ter than one percent (l3| . It is therefore important to 
see if an analogous level of precision in the prediction of 



the Casimir force can be obtained at the theoretical level. 
A precise determination of the theoretical error is indeed 
as important as reducing the experimental error, in order 
to address controversial questions that have emerged in 
the recent literature on dispersion forces, regarding the 
influence of free charges on the thermal correction to the 
Casimir force [l^ . 

Addressing the theoretical error in the magnitude of 
the Casimir force is indeed difficult, because many phys- 
ical effect must be accounted for. However, it has re- 
cently been pointed out [l9j that perhaps the largest the- 
oretical uncertainty results from incomplete knowledge 
of the optical data for the surfaces involved in the ex- 
periments. On one hand, the large variability depending 
on the preparation procedure, of the optical properties 
of gold coatings, routinely used in Casimir experiments, 
makes it necessary to accurately characterize the coatings 
actually used in any experiment. On the other hand, 
even when this characterization is done, another prob- 
lem arises, because for evaluating the Casimir force one 
needs to determine the electric permittivity e(i£) of the 
coatings for certain imaginary frequencies i£. This quan- 
tity is not directly accessible to any optical measurement, 
and the only way to determine it is via exploiting disper- 
sion relations, that permit to express e(i£) in terms of 
the measurable values of the permittivity e(w) for real 
frequencies lu. When doing this, one is faced with the 
difficulty that optical data are necessarily known only in 
a finite interval of frequencies CL> m ; n < to < to max . This 
practical limitation constitutes a severe problem in the 
experimentally relevant case of good conductors, because 
of their large conductivity at low frequencies. With the 
standard forms of dispersion relations Eq. ^ and Eq. 
(|9]), one finds that for practical values of uj m - m and w max , 
low frequencies less than w m i n and/or large frequencies 
larger than w max give a very large contribution to e(i£). 
In order to estimate e(i£) accurately, one is then forced 
to extrapolate available optical data outside the experi- 
mental region, on the basis of some theoretical model for 
e(cj). Of course, this introduces a further element of un- 
certainty in the obtained values of e(i£), and the resulting 
theoretical error is difficult to estimate quantitatively. 

In this paper we have shown that this problem can be 
resolved by suitably modifying the standard dispersion 
relation used to compute e(i£), in terms of appropriate 
analytic window functions f(z) that suppress the contri- 
butions both of low and large frequencies. In this way, 
it becomes possible to accurately estimate e(i£) solely on 
the basis of the available optical data, rendering unneces- 
sary any uncontrollable extrapolation of data. We have 
checked numerically the performance of simple choices 
of window functions, by making a numerical simulation 
based on an analytic fit of the optical properties of gold, 
that has been used in recent experiments on the Casimir 
effect [l2j]. We found that already very simple forms 
of the window functions permit to estimate the Casimir 
pressure with an accuracy better than one part in a thou- 
sand, on the basis of reasonable intervals of frequencies 
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for the optical data. It would be interesting to apply 
these methods to the accurate optical data for thin gold 
films quoted in Ref. 

Before closing the paper, we should note that the rele- 
vance of the sample-to-sample dependence of the optical 
data observed in [r| for the theory of the Casimir effect 
has been questioned by the authors of Ref. [l2|, who 
observed that this dependence mostly originates from 
relaxation processes of free conduction electrons at in- 
frared and optical frequencies, due for example to differ- 
ent grain sizes in thin films. The main consequence of 
these sample-dependent features is the large variability 
of the Drude parameters, extracted from fits of the low- 
frequency optical data of the films, which constitutes the 
basic source of variation of the computed Casimir force 
reported in Ref. [l9|. According to the authors of Ref. 
[12] , relaxation properties of conduction electrons in thin 
films, described by the fitted values of the Drude pa- 
rameters, are not relevant for the Casimir effect. Indeed, 
according to these authors the quantity e(u) to be used in 
Lifshitz formula should not be understood as the actual 
electric permittivity of the plate, as derived from opti- 
cal measurements on the sample, but it should be rather 
regarded as a phenomenological quantity connected to 
but not identical to the optical electric permittivity of 
the film. The ansatz offered by them for e(w) is dubbed 
as generalized plasma model, and following Ref. [HJ we 



prescribes that only the core-electron contribution e c (ui) 
should be extracted from optical data of the film. On 
the contrary, and more importantly, according to Ref. 
Il2j the value of the plasma frequency uj p to be used in 
Eq. (f3Tj) should be the one pertaining to a perfect crys- 
tal of the bulk material, and not the one obtained by a 
Drude-model fit of the low-frequency optical data of the 
film actually used in the experiment. The justification 
provided for this choice of the plasma frequency by the 
authors of Ref. [l2j is that the contribution of conduc- 
tion electrons to the Casimir force should depend only 
on properties determined by the structure of the crys- 
tal cell, which are independent of the sample-to-sample 
variability determined by the peculiar grain structure of 
the film, reported in Ref. [l!|. It should be noted that 
for gold, the value of the plasma frequency advocated in 
[l^ . uip = 9 eV/h, is much higher than the fit values 
quoted in Ref. 19], which range from 6.8 to 8.4 eV/Ti. 
As a result, the approach advocated in Ref. [l2| leads to 
larger magnitudes of the Casimir force, as compared to 
the values derived in Ref. [l9j], with differences ranging, 
depending on the sample, from 5 % to 14 % at 100 nm. 
There is no room here to further discuss the merits and 
faults of these approaches, and we refer the reader to 12 1 
for a thorough analysis. It is fair to note though that a 
series of recent experiments by one experimental group 
[IsH appears to favor the generalized plasma approach, 



denote it as e gp (w). This quantity is a semianalytical and to rule out the more conventional approach based on 



mathematical construct, defined by the formula: 
£ g P M = e c {uj) ~uj 2 p /uj 2 , 



actual optical data followed in Refs. 



approa 



(31) 



where t c {uj) represents the contribution of core electrons, 
while the term proportional to the square of the plasma 
frequency uj p describes conduction electrons. The most 
striking qualitative feature of this expression is the ne- 
glect of ohmic dissipation in the contribution from con- 
duction electrons, but this is not all. Indeed, the ansatz 



The future will tell what is the correct description. 
In the meanwhile, we remark that whatever approach is 
followed, the methods proposed in this paper may prove 
useful to obtain more reliable estimates of the Casimir 
force for future experiments. 
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